Compressibility and Electronic Structure of MgB2 up to 8 GPa 



T. Vogt and G. Schneider 
Physics Department, Brookhaven National Laboratory, Upton, NY 11973-5000, USA 



o 
o 

(N 
D 

(N 



C/3 



J. A. Hriljac 

School of Chemistry, The University of Birmingham, Edgbaston, Birmingham B152TT, United Kingdom 

G. Yang and J.S. Abell 

School of Metallurgy and Materials, The University of Birmingham, Edgbaston, Birmingham B152TT, United Kingdom 

(February 1, 2008) 

The lattice parameters of MgB2 up to pressures of 8 GPa were determined using high-resolution 
x-ray powder diffraction in a diamond anvil cell. The bulk modulus, Bq, was determined to be 
151±5 GPa. Both experimental and first-principles calculations indicate nearly isotropic mechanical 
behavior under pressure. This small anisotropy is in contrast to the 2 dimensional nature of the boron 
p states. The pressure dependence of the density of states at the Fermi level and a reasonable value 
for the average phonon frequency {uj) account within the context of BCS theory for the reduction 
of Tc under pressure. 
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The recent discovery of superconductivity in MgB2 
close to 40K [Q has pushed Tc beyond what was thought 
to be possible in intermetaUic compounds within the 
context of BCS theory [^jsj. However, the boron iso- 
tope effect was measured to be consistent with MgB2 be- 
ing a conventional phonon-mediated BCS superconduc- 
tor Uj.The temperature-dependent electrical resistivity 
in the normal state was observed to obey a power 
law by Finnermore, et al. whereas Jung et al. ||] 
found a behavior to be more appropriate. Within 
the context of hole superconductivity [Q it was proposed 
that a decreasing B-B intraplane distance should increase 
Tc- However, ac magnetic susceptibility measurements 
revealed that actually decreases jsj. The response of 
the crystal and electronic structures of MgB2 to pres- 
sure is not only important to distinguish a conventional 
from a non-BCS mechanism, but it might also give valu- 
able clues for guiding chemical substitutions. Slusky, et 
al. 1^ have established the loss of superconductivity in 
Mgi_2;Ala;B2 and conjecture that MgB2 is on the edge of 
a structural instability at slightly higher electron concen- 
trations. These authors furthermore indicate the presence 
of a two phase region between 0.1 0.25. To inves- 

tigate the intrinsic electronic properties of MgB2 and to 
avoid complications from possible extrinsic effects such 
as compositional fluctuations and phase separations in 
substituted systems, we have determined the bulk com- 
pressibility and compared it to the one derived from first- 
principles calculations based on Density Functional The- 
ory (DFT). 

The sample, obtained from Alfa, was of excellent qual- 
ity without any impurities visible in a high-resolution 
x-ray powder diffraction pattern taken at room temper- 



ature. All reflections could be indexed to the AlB2-type 
structure in spacegroup P6/mmm with a=3. 08589(1) A 
and 0=3.52121(3) A at room temperature and ambient 
pressure. The magnetic susceptibility showed a sharp 
transition at 38. 2K (Figure 1). The compressibility up 
to 8 GPa was measured using a Merrill-Bassett diamond 
anvil cell. The lattice parameters (Table 1) were de- 
termined using the high-resolution powder diffractome- 
ter X7A at the National Synchrotron Light Source at 
Brookhaven National Laboratory. An asymmetrically 
cut triangular shaped Si (220) monochromator was cylin- 
drically bent to focus a photon beam with a wave- 
length of 0.68613(2) A down to 200 ^m. A steel gas- 
ket was pre-indented and a 200 /im hole drilled in the 
center. The sample and several small ruby chips were 
loaded along with the pressure transmission fluid of 
methanol:ethanol:water (ratios 16:3:1). Pressure was de- 
termined by measuring the shift of the Rl emission line 
of ruby using an argon ion laser pO[ |. Diffraction pat- 
terns up to 29 = 35° were recorded using a position- 
sensitive detector Further experimental details are 
described in ||l^ . There were no indications of any split- 
ting or significant broadening of the Bragg refiections 
indicative of a possible phase transition. The width of 
the Rl emission line of ruby at higher pressures also in- 
dicated that the sample was well within a hydrostatic 
pressure regime. The unit cell volumes, V , normalized 
to the one at ambient pressure, Vo, at various pressures 
P were fitted to a first-order Murnaghan equation of state 
V = Vo(l + B'gP/Boy^^^"' where Bq is the bulk modu- 
lus at ambient conditions and Bq/ = 4. We obtained a 
value of 151 ± 5 GPa. For comparison, Bq of YNi2B2C, 
a related borocarbide, was measured to be 200 GPa [n3|. 
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We calculated the electronic structure and mechani- 
cal properties within DFT using the full-potential 
linearized augmented plane wave method FLAPW p5| . 
For the exchange-correlation potential we used the Gen- 
eralized Gradient Approximation of Perdew, et al. [ p^ . 
A relatively large basis set was required to accurately 
calculate small changes in the electronic structure and 
lattice parameters with applied pressure 10 . Calculated 
equilibrium parameters are a — 3.089A, c = 3.548A and 
c/a = 1.149. They differ from the experimental values by 
0.1%, 0.8% and 0.7% respectively. The calculated bulk 
modulus of Bq = 139 ±10 GPa is in good agreement with 
the measured value of 151 ± 5 GPa. Both calculated and 
experimental lattice parameters as a function of volume 
are shown in Fig. 2. The agreement between theory and 
experiment is excellent. Over the range of pressure con- 
sidered, the c/a ratio is essentially constant, indicative 
of the 3 dimensional character of MgB2. The small ob- 
served anisotropy of the compressibility as given by the 
observed reduction in the c/a ratio with decreasing vol- 
ume is reproduced by our calculations, although the rate 
of change is overestimated. 

The superconducting properties of MgB2 as a function 
of pressure have been investigated by Lorenz, et al. 
who found that the transition temperature Tc decreases 
linearly and reversibly at a rate of —1.6 K/GPa. Within 
BCS theory the pressure dependence of Tc can be calcu- 
lated from the the McMillan formula 



Tc = (w) exp 



-(1 + A) 



X~ fi* - Xfi* 



where (uj) is an average phonon frequency, /i* denotes the 
Coulomb pseudopotential, A oc A^(0)/ (tj^) is the electron- 
phonon interaction parameter and A^(0) denotes the den- 
sity of states at the Fermi energy . 

The calculated change in the density of states of MgB2 
for a pressure of w 8 GPa is shown in Figure 3. We ob- 
serve the expected increase in bandwidth, but the overall 
changes are very small and in particular A^(0) decreases 
by only 3% from 0.70 states/eV to 0.68 states/eV. As- 
suming that Tc continues to decrease linearly at higher 
pressures we can estimate that at 8 GPa, Tc « 26K which 
is 30% lower than Tc at ambient pressure. Comparing 
this large decrease in Tc to the minute change in A^(0), it 
becomes apparent that the pressure dependence of A^(0) 
alone cannot account for the large rate with which Tc 
decreases with applied pressure. Taking our value for 
d\iiN{0)/dP « -0.004/GPaand assuming n* = 0.1 and 
A = 1.0, which is the average of published calculations 
and estimates ranging from 0.7 to 1.4 jl^Jl^, we find that 
in order to reproduce the experimentally observed pres- 
sure dependence of Tc, dlTi{uj'^)^/'^ /dP w -1-0.014/ GPa is 
needed. This value appears reasonable. 

The volume dependence of Tc within BCS theory has 
been characterized using 



d\n{Tc/eD) , Qd d\nX , 9^ 

= in = In G9 

dlnV Tc d\nV Tc 

where ip is a material dependent parameter describing 
the volume dependence of the electron-phonon coupling 
parameter and Od is the Debye temperature. For su- 
perconducting sp metals one generally finds « 2.5 [po[ . 
Using for the Debye temperature the recently determined 
value of = 800 K 111, we find (p = 2.4 for MgBa, a 
value not inconsistent with BCS theory. 

The importance of the mostly 2-dimensional B p^y 
states for the superconductivity of MgB2 has been 
pointed out by several authors |l^Jl^. The change with 
pressure in the partial density of states (PDOS) of B p 
states separated into Pxy and Pz character is shown in 
Fig. 4. The B pz PDOS remain essentially unchanged, 
but the B Pxy PDOS shows an overall shift to lower en- 
ergies, leading to a loss of states in the energy range 
from just below the Fermi energy to « —2 eV. In gen- 
eral one expects the PDOS inside a fixed sphere to in- 
crease with pressure. The qualitatively different behavior 
of the B Pxy PDOS in MgB2 can be attributed to the 
2-dimensional character of these B p^y states. However, 
the B Pxy PDOS at the Fermi energy does not change 
significantly. 

In conclusion, we have experimentally and theoreti- 
cally determined the bulk modulus of MgB2 and lattice 
parameters up to pressures of 8 GPa. Despite an overall 
isotropic behavior of MgB2, there are small anisotropics 
in the mechanical properties. The electronic structure 
under pressure reveals the 2-dimensional character of the 
B Pxy states whose occupancy is altered by pressure in 
contrast to the B pz states. 

Furthermore, the pressure dependence of the density of 
states at the Fermi level and a reasonable value for the 
average phonon frequency (w) account for the reduction 
of Tc under pressure within the context of BCS theory. 

We find that the material dependent constant ip de- 
scribing the volume dependence of the electron-phonon 
coupling parameter agrees well with BCS theory for an 
sp metal. 

We gratefully acknowledge discussions with Mike 
Weinert and Myron Strongin. We also would like to 
thank J. Hu from the Geophysical Laboratory, Carnegie 
Institution of Washington for being able to use the laser 
system for pressure determination. The work was sup- 
ported by the Division of Materials Science, US De- 
partment of Energy, under Contract NO. DE_AC02- 
98CH10886 and a grant of computer time at the National 
Energy Research Scientific Computing Center (NERSC). 



2 



[1] 

[2] 
[3] 



[4] 

[5] 

[6] 

[7] 
J8L 



[9] 



[10] 

[11] 
[12] 

[13] 

[14] 



[15] 



[16] 



[17] 



[18] 
[19] 
[20] 



[21 



J. Akimitsu, Symposium on Transition Metal Oxides, 
(Sendai, Japan, January 10, 2001). 
W.L. McMillan, Phys. Rev. 167, 331 (1968). 
M.L. Cohen, P. W. Anderson in Superconductivity in d- 
and f-band metals, edited by D.H. Douglass, AIP Con- 
ference Proceedings No 4, (AIP, New York 1972), page 
17. 

S. L. Bud'ko, G. Lapertot, C. Petrovic, C. E. Cunning- 
ham, N. Anderson, and P. C. Canfield PRL 86,1877 
(2001). 

D.K. Finnermore, J.E. Ostenson, S.L. Bu d'ko, G. Laper- 
tot and P.C. Canfield [eond-mat 0102114 
C. U. Jung, M. S. Park, W. N. Kang . M. S. Kim, S. Y. 
Lee, and S. L Lee |:ond-mat 01022"lj r 



J. E. Hersch sond-mat 0102115 

B. Lorenz, R. L. Meng, and C. W. Chu, cond- 



mat/0102264. 



Regan, M. A. Hayward, 
S. Loureiro, M. K. Haa s, 
cond- /010226^ . 



J. S. Slusky, N. Rodago, K. A. 
P. Khalifah, T. He, K. Inumaru, 

H. W. Zandbergen, and R. J. Cava 
G. J. Permarini, S. Block, J. D. Barnett, and R. A. For- 
man J.Appl. Phys. 46, 2774 (1975). 

G. Smith Synchr. Rad. News 4, 24 (1991). 
A. P. Jephcoat, L.W. Finger, and D.E. Cox, High Pres- 
sure Research 8, 667 (1992). 

S. Meenakshi, et al., Physica B 223-224, 93 (1996). 

P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 

(1964) , W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 

(1965) . 

E. Wimmer, H. Krakauer, M. Weinert, and A. J. Free- 
man, Phys. Rev. B 24, 864 (1981); H. J. F. Jansen and 
A. J. Freeman, Phys. Rev. B 30, 561 (1984). 
J. P. Perdew, K. Burke, and M. Enzerhof, Phys. Rev. 
Lett. 77, 3865 (1996); J. P. Perdew, K. Burke, and M. 
Enzerhof, Phys. Rev. Lett. 78, 1396(E) (1997); J. P. 
Perdew, K. Burke, and M. Enzerhof, Phys. Rev. Lett. 
80, 891 (1998); Y. Zhang and W. Yang, Phys. Rev. Lett. 
80, 890 (1998). 

The parameters used in all calculations are Rmt ~ 

I. 6(2.7) a.u. for B (Mg) atoms. Potential and charge den- 
sity were expanded up to Imax = 8 and Gmax = 18 a.u.~^ . 
The wave-function cutoff was chosen to be RMrkmax ~ 
8.8 leading to ~ 180 basis functions per atom. 
J. 



:ond-mat/0101446 



Kortus et al. 
J. M. An and W. E. Pickett, |cond-mat/0102391 . 
M. Levy, J. L. Olsen in Physics of High Pressures and 
the Condensed Phase edited by A. van Itterbeek, North 
Holland, 1965, Ch. 13. 

R. K. Krem er, B. J. Gibson, and K. Ahn cond 



mat/0102432 



TABLE L Measured lattice parameters for MgB2 as a func- 
tion of applied pressure. The error in pressure is estimated to 
0.2 GPa. 



pressure (GPa) volume (A=* 



(A) 



(A) 



c/a 



0.0 


29.0391 


1.17 


28.8494 


2.14 


28.5835 


3.05 


28.5017 


4.07 


28.2994 


5.09 


28.0519 


6.53 


27.8583 


8.02 


27.8228 



3.08589(3) 
3.08017(28) 
3.07149(21) 
3.06709(49) 
3.06349(15) 
3.05449(19) 
3.04972(31) 
3.04843(14) 



3.52121(9) 
3.51121(29) 
3.49854(52) 
3.48857(58) 
3.48188(25) 
3.47180(24) 
3.45863(24) 
3.45715(28) 



1.1411 
1.1399 
1.1390 
1.1374 
1.1366 
1.1366 
1.1341 
1.1341 
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FIG. 1. ac susceptibility of MgB2 as a function of temper- 
ature. 
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FIG. 2. c/a ratio as a function of unit-cell volume from 
both experiment and calculation. 
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FIG. 4. Single Boron p Partial Density of States (PDOS) 
(top panel) for two sets of lattice parameters (see Figure 3), 
corresponding to Gpa (full lines) and « 8 GPa (dashed 
lines). The change in the B p^y PDOS (full line) andp^ PDOS 
(dashed line) is shown in the bottom panel. The energy zero 
is the Fermi energy. 
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FIG. 3. Total and Partial Densities of States (PDOS) of 
MgB2 for two sets of lattice parameters: ( a = 3.09 A, 
c = 3.55 A) (full lines) and ( a = 3.05 A, c = 3.46 A) (dashed 
lines). The difference in lattice parameters corresponds to 
a pressure of ~ 8 GPa. The energy zero is the Fermi en- 
ergy. The PDOS is calculated inside the Muffin tin spheres ( 
Rmt = 1.6 a.u. and 2.7 a.u. for B and Mg ). 
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